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The interpretation of seismoelectric observations involves the dynamic electrokinetic cou- 
pling, which is related to the streaming potential coefficient. We describe the different 
models of the frequency-dependent streaming potential, mainly the Packard's and the 
Pride's model. We compare the transition frequency separating low-frequency viscous 
flow and high-frequency inertial flow, for dynamic permeability and dynamic streaming 
potential. We show that the transition frequency, on a various collection of samples for 
which both formation factor and permeability are measured, is predicted to depend on 
the permeability as inversely proportional to the permeability. We review the experi- 
mental setups built to be able to perform dynamic measurements. And we present some 
measurements and calculations of the dynamic streaming potential. 
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1 INTRODUCTION 

Electrokinetics arise from the interaction between the rock matrix and the pore water. There- 
fore electrokinetic phenomena are often observed in aquifers, volcanoes, and hydrocarbon 
or hydrothermal reservoirs. Observations show that seismoelectroma gnetic signals associ 



ated to earthquakes can be induced by electromagn etic induction 



Matsushima et al.l 120021 ) or by electrokinetic effect (ITakeuchi et al. 



(Ho nkura et al. 



1998 



2009: 



Fenoglio et al. 



19951 ). The electrokinetic phenomena are due to pore pressure gradients leading to fluid 
flow in the porous media or fractures, and inducing electrical fields. These electrokinetic 
effects are associated to the electrical double layer which was originally described by Stern. 
The electrokinetic signals can be induced by global displacements of the reservoir fluids 
(streaming potential) or by the propagation of seismic waves (seismoelectromagnetic effect). 
As soon as these pressure gradients have a transient signature, the dynamic part of the 
electrokinetic coupling has to be taken into account by introducing the dependence on fluid 
transport properties. 

It is generally admitted that two kinds of seismoelectromagnetic effects can be ob- 
served. The dominant contribution, commonly called "coseismic", is generated close to 
the receivers du ring the passage o f seismic waves. The second kind, so called "interfa- 



cial conversion" (IDupuis et al. 



20091 ). is very similar to dipole radiation and is generated 
at physico-chemical interfaces due to strong electrokinetic coupling discontinuities. This in- 
terface conversion is often perceived to have the potential to detect fine fluids transitions 
with higher resolution than seismic investigations, but in practice, signals are often masked 
by electromagnetic disturbances, especially when generated at great depth. 

Nevertheless recent field studies have focused on the seismo-electric conversions linked 



to electrokinetics in order to investigate oil and gas reservoirs fjThompson 



drau 



2007 



i c reservoirs (IDupuis fc Butler 



Haines et al. 



2007a 



2006 



Dupuis et al. 



Strahser et al. 



2011 



2007; 



et al. 



Du puis et a 



Garambois fc Dietrich 



2009 



200 51) or hy- 



Strahser et al 



200 lh . It has been 



shown using these investigations that not only the depth of the reservoir can be deduced, 
but also the geometry of the reservoir can be imaged using the amplitudes of the electro- 



seismic signals ( iThompson et al. 
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20071 ). Moreover fractured z ones can be detected and per 



mcabi 



2005 



ity can be measured using; seismo-elect rics in borehole ( iSinger et al. 



2005 



Pain et al. 



Mikhailov et al. 



2000 



Jouniaux 



201ll ). This method is especially appealing to hy- 



drogeophysics for the detection of subsurface interfaces induced by con trasts in permeabil 



ity, in porosity, or in electrical properties (salinity and water content) ( ISchakel et al 



2011 



Schakel fc Smeulders 



2010 



Garambois fc Dietrich 



20021 ). 



The a nalytical inte rpretation of the seismoelec tromagnetic phenomenon has been de 



scribed by 



Pridd (Il994t ). by connecting the theory of 



Biot 



(119561 ) for the seismic wave propa- 



gation in a two phases medium with Maxwell's equations, using dynamic electrokinetic cou- 



plings. The seism oelectromagnetic conversions have been mode 



saturated media ( 



2002 



Gao k Hu 



Haartsen fc Pride 



1997 



Haartsen et al. 



1998 



ed in homogeneous or layerec 



Garambois Sz Die 



20101 ) with applications to reservoir geophysics (jSaunders et al 



rich 



20061 ) 



2001 



Theoretical developments showed that the electrical field induc ed by the P- waves prop- 



agation is related to the acceleration ( IGarambois fc Dietrich 



20011). The electrokinetic cou- 



pling is created at the interface between grains and water, when there is a relative mo- 
tion of electrolyte ions with respect to the mineral surface. Thus, seismic wave propaga- 
tion in fluid-filled porous media generates conversions from seismic to electromagnetic en- 
ergy which can be observed at the macroscopic scale, due to this electrokinetic coupling 
at the pore scale. The seismoelectric coupling is directly dependent on the fluid conduc- 



tivity, the fluid density and the electric dou ble-layer (the electrical inter: 



grains and the water) (see the tutorial by ( jJouniaux &: Ishido 



this issue 



ace between the 
), in this special 



issue "Electrokinetics in E arth Science s ' for m or e details). For more details on the surface 



complexation reactions see 



Davis et al. 



f|l978f ) or 



Guichet et al 



quantified in the broad band by a dynamic coupling (IPride 



(120061) . It can be accurately 



19941 ) which can be linked 



in the low frequency l imit to the steady-state 
ied in p orous media (|lshido fc Mizutani 



1995a 



b. 



1997 



Jouniaux et al 



1994, 



1999 



198 



2000 



s treaming potentia l coefficient largely stud- 



Pozzi & Jouniaux 



Guichet et al. 



1994 



2003 



2006 



Jouniaux fc Pozzi 



Jaafar et al 



2009; 



Jouniaux et al. 



2009 



Vinogradov et al 



20101 ; 



Jackson 



2010 



Allegre et al. 



20101 ). 
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seismoelectric conversions ( 


Vligunov & Kokorev 


1977: 


Chandler 


1981 




A 


ironov et al. 


1994; 


Jiang et al. 


1998: 


Zhu 


et al. 


1999. 


2000 


; Zhu 


& Toksoz 


2003; 


Chen & Mu 


2005: 


Bordes et al. 


2006; 


Block & Harris 


20o4 


Zhu et al. 


2008; 


Bordes et al. 


2008 


). These papers describe the 



laboratory studies performed to investigate this dynamic coupling. An oscillating pore pes- 
sure must be applied to a rock sample, and because of the relative motion between the rock 
and the fluid, an induced streaming potential can be measured. Depending on the oscillating 
frequency of the fluid, the fluid makes a transition from viscous dominated flow to inertial 
dominated flow. As the frequency increases, the motion of the fluid within the rock is de- 
layed and larger pressure is needed. In order to know the dynamic coupling, both real and 
imaginary part of the streaming potential must be measured. 

2 FROM DYNAMIC STREAMING POTENTIAL TO 
SEISMOELECTRO MAGNETIC COUPLING 

The steady-state streaming potential coefficient is defined as the ratio of the streaming 
potential to the driving pore pressure: 



AP rjCTf 

which is called the Helmholtz-Smoluchowski equation, where (jf, e and rj are the fluid con- 
ductivity, the dielect ric constant of the fluid, and the fluid dynamic viscosity respectively 



(see the tutorial by (Uouniaux fc Ishido 



this issue! ) ). In this formula the surface electrical 



conductivity is neglected compared to the fluid electrical conductivity. The potential ( is 
the electrical potential within the double-layer on the slipping plane. Although the zeta 
potential can hardly be modeled for a rock and although it can not be direclty measured 
within a rock, the steady-state streaming potential coefficient can be measured in labora- 
tory, by applying a fluid pressure difference ( AP) and by measuring the induced streaming 



electric potential (AV) ( 1 Jouniaux et al. 



2000 



Guichet et al. 



2003 



2006 



Allegre et al. 



2010 



Frequency- dependent streaming potentials: a review 5 



201ll ). The el ectrical potential C itself depends on fluid composition and pH, and the water 



conductivity ( 



Davis et al. 



1978 



Ishido fc Mizutanil 



1981 



Lome et al. 



1999 



Jouniaux et al. 



201oL 



2000 



Guichet et al 



2006 



Jaafar et al 



2009 



Vinogradov et al 



2010 



Allegre et al 



2.1 Packard's model 



Packard! (119531 ) proposed a model for the frequency-dependent streaming potential coefficient 
for capillary tubes, assuming that the Debye length is negligible compared to the capillary 
radius, based on the Navier-Stokes equation: 

AV(u) 



C s0 (uj) 



APM 



)(- 



a 



lUlpf 



(2) 



v vv v v 

where u is the angular frequency, a is the capillary radius, J% and Jo are the Bessel 
functions of the first order and the zeroth order, respectively, and pf is the fluid density. 
The transition angular frequency for a capillary is: 



More recently 



i r 



Pf a 



(3) 



Reppert et al.l ( 1200 ll ) used the low- and high-frequency approximations 



of the Bessel functions to propose the following formula, which corresponds to their eq.26 
corrected with the right exponents —2 and —1/2: 



C s nM 



Wf 



with the transition angular frequency 



UJ r 



~2_ plfl L, 

87] 



Pft 



-2 



(4) 



(5) 



and s howed that this model was not very different from the model proposed by 



Packard 



dl953j). 

The complete development relating the Biot's theory and the Maxwell's equations has 
been published by Pride in 1994. 
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2.2 Pride's model 



Pridd ( 11994] ) derived the equations governing the coupling between seismic and electromag- 
netic wave propagation in a fluid-saturated porous medium from first principles for porous 
media. The following transp ort equation s express the coupling between the mechanical and 



electromagnetic wavefields [( jPride 



19941 ) equations (174), (176), and (177)]: 



cr(u)E + L(u) (-Vp + iu 2 p f u s ) 



(6) 



— iuw = L(u/)E 



V 



(-Vp + iuj 2 p f u s ) 



(7) 



In the first equation the macroscopic electrical current density J is the sum of the average 
conduction and streaming current densities. The fluid flux w of the second equation is 
separated into electrically and mechanically induced contributions. The electrical fields and 
mechanical forces that create the current density J and fluid flux w are, respectively, E 
and (— Vp + iu 2 pfU s ), where p is the pore-fluid pressure, u s is the solid displacement, and 
E is the electric field. The complex and frequency-dependent electrokinetic co upling L(u) 



which describes the coupling between the seismic and electromagnetic fields (jPride 



Reppert et al. 



1994 



20011 ) is the most important parameter in these equations. The other two 
coefficients, cr(oj) and k(oj), are the electric conductivity and dynamic permeability of the 
porous material, respectively. 

The seismoelectric coupling that describes the co upling betwe en the seismic and electro- 



magnetic fields is complex and frequency-dependent 



Pride! (jl994f ): 



L(u) = Lr 



1 -^?( 1 - 2 l) 2 ( 1 -' 3/2d \/f) : 



(8) 



where L Q is the low frequency elect rokinetic coup 



a porous- m ateria 
(detailed in 



geome try term (IJohnson et al 



ing, d is related to the Debye-length, A is 



Pridd (11994( 1) 



19871 ). and m is a dimensionless number 



The transition angular frequency u c separating low-frequency viscous flow and high- 
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frequency inertial flow is defined as: 

<* = (9) 
where is the porosity, k Q is the intrinsic permeability, is the tortuosity. 

2.3 Further considerations 

The low-frequency electrokinetic coupling L Q is related to the steady-state streaming poten- 
tial coefficient C s q by: 

L = -C s0 a r (fO) 

where o~ r is the rock conductivity The electrokinetic coupling L(oS) can be estimated by 
considering that steady-state models of C s0 can be applied to the calculation of L . When 
writting o r = o~f/F with surface conductivity neglected, the steady-state electrokinetic 
coupling can be written as: 

* = (») 

We can see that the steady-state electrokinetic coulping is inversely proportional to the 
formation factor. 

The transition angular frequency separating viscous and inertial flows in porous medium 
can be rewritten by inserting = cf) F with F the formation factor that can be deduced 
from resistivity measurements using Archie's law, as: 

^■44 (12) 

where F is the formation factor that can be deduced from resistivity measurements using 
Archie's law. 

Since the per meability and t he formation factor are not independent, but can be related 



by k = CR 2 /F flPaterson 



19831 ) with C a geometrical constant usually in the range 0.3-0.5 
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and R the hydraulic radius, the transition angular frequency can be written as: 

(13) 

The equation [13] shows that the transition angular frequency in porous medium is inversely 
proportional to the square of the hy draulic radius. 



Recently IWalker fc Glover! (|2010[ ) proposed a simplified equation of Pride's development 
assuming that the Debye length is negligible compared to the characteristic pore size, and 
assuming the parameter: 



leading to the equation: 



m 



Liu) 



A 



eff 



2i- 



A 



eff 



(14) 



(15) 



with r e ff the effective pore radius, and a transition angular frequency 



577 



Pf r lff 



(16) 



Garambois fc Dietrichl ( 1200 ll ) studied the low frequency assumption valid at seismic fre- 
quencies, meaning at frequencies lower than the Biot's frequency separating viscous and 
inertial flows and gave the coseismic transfer function for low frequency longitudinal plane 
waves. In this case, and assuming the Biot's moduli C « H, they showed that the seismo- 
electric field E is proportional to the grain acceleration: 



E 



L .. < 

~ PfU = 

a r 7/07 



p f u 



(17) 



Equations [T71 [10] and [T] show that transient seismo-electric magnitudes will be affected 
by the bulk density of the fluid, and the streaming potential coefficient which is inversely 
proportional to the water conductivity and proportional to the zeta potential (which depends 
on the water pH). 
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2.4 The electrokinetic transition frequency compared to the hydraulic's one 

he theory of dynamic permeability in porous media has been studied by many authors 



(lAuriault et al.l 



1222). 



1985 



IJohnson et al. 



1987 



Sheng k Zhoulll988 



Sheng et al. 



1988 



Smeulders et al 



The frequency behavior of the permeability is given by 



Pride! (11994! ) by: 



kjuj) 
k 



u 4 



n -1 



(jj r 



(18) 



The transition angular frequency for a porous medium is the same as eq. [91 



Charlaix et al. 



(119881 ) measured the behavior of permeability with frequency on capillary tube, glass beads 



and crushed glass. The dynamic permeability is constant up to the transition frequency 
above which it decreases, and the more permeable the sample is, the lower the transition 



frequency is. Ot 



ler m easurements have been performed on glass beads and sand grains 



(jSmeulders et al 



1992). The transition frequency (/ c = u c /2n) varies from 4.8 Hz to 149 
Hz for samples having permeability in the range 10 -8 to 10~ 10 m 2 (see Tabled]), which are 
extremely high permeabilities. 

The transition frequency indicates the beginning of the transition for both the perme- 
ability and the electrokinetic coupling. However the transition behavior and the cuttoff 
frequency are different between permeability and electrokinetic coupling (eq. Eland eq JTBl . 
both depending on the pore-space geometry term m but in different manner. 

We calculated the predicted transition frequency f c = u c /2tt with u c from eq. fT2l with 
7] = 10~ 3 Pa.s and pf = 10 3 kg/m 3 . T he other parameters F and ko are measured from 
different authors cited in lBernaba (ll99ll ) (see Table [2]). We also calculated the parameters for 
four F ontainebleau sandst one samples. It has been shown for these samples that F = 0~ 2 01 



(from iRuffet et al. 



(Il99ll )) and that ko = a(j) n with different values for n according to the 



porosity. The following laws were cho sen: kn = 1.66x 



for cf) ranging between 8 and 25% (IBourbie et al. 



CT 



for 6 < 6% and k n = 2.5x10 



-HU3 



19871 ). We can see that the transition 



frequencies are of the order of kHz and MHz and no more from 0.2 to 150 Hz as measured or 
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calculated on glass beads, sand grains, crushed glass or capillaries. We plotted the results of 
the transition frequency as a function of the permeability on these various samples in Fig. 
[TJ Although the formation factor is not constant with the permeability, it is clear that the 
transition frequency is inversely poportional to the permeability as: 



log w (f c ) = -0.78log 10 (k) - 5.5 (19) 

and varies from about 100 MHz for 10~ 17 m 2 to about 10 Hz for 10 -8 m 2 , so by seven 
orders of magnitude for nine orders of magnitude in permeability. 



3 EXPERIMENTAL APPARATUS AND PROCEDURE 



Several experimental setups were proposed to provide the sinusoidal pressure variations. 
The first experimental apparatus proposed a sinusoidal motion delive red by a sylpho n 



bellows which was driven by a geophone-type push-pull driver (Fig. 121 from iPackardl (119531 )) 
The low frequency oscillator (0.01 Hz to 1 kHz) was used for opera tion of the pu sh-pull 



geophone driver. Similar setups were proposed by 



Thurston! (1952b) 



1955) , so tha t frequency of t his kind o 



([Packard 



More recently 



1953) and 10-700Hz f 



Schoemaker et al 



rurston 



source was 1-400 Hz rtCooke 



iff. [3) and 



Cooke 



1955|), 20-200 Hz 



1952bl ). The induced pressure was up to 2 kPa. 



(120071 ) used a so-called Dynamic Darcy Cell (DCC) with 
a mechanical shaker connected to a rubber membrane leading to a frequency range for the 
oscillating pressure 5 to 200 Hz. The sinusoidal fluid flow was also applied b y a displacement 



Groves fc Sears 



piston pump directly connect ed to the electrodes chambers (fig. H] from 
(119751 ); ISears fc Groves! (119781 )). The piston was mounted on a S cotch Yoke drive attached 



to a controllable speed AC motor (jCerda fc Non-Chhom 



1989|). The 



frequency range o f 



Pengra et al. 



(1999) 



this source was then 0.4Hz to 21 Hz and the pressure up to 15 kPa. 
used a piston rod attached to a loudspeaker driven by an audio power amplifier (Fig. [5]). 
They performed measurem ents up to 100Hz, wi th an applied pressure of 5 kPa RMS. More 



recently it was proposed by 



Reppert et al. 



(]200ll ) to use an electromechanical transducer (fig. 
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El), and these authors cover ed a frequency range 1-500 Hz. The vi brating excit e r pro posed 



by 



Schoemaker et al. 



( 120081 ) was used from 5Hz to 200Hz. Recently iTardif et al. 



(120111 ) used 



an electromagnetic shaker operating in the range 1Hz to 1kHz anc 

requencies have 



up to 200Hz. Higher 



2005 



Block fc Harris 



2006 



Zhu et al. 



3een investigated (IZhu et al. 



proy i ded measurements 



1999 



2000 



Chen k Mu 



20081 ) for the detection of the interfacial conversions. 



The electromagnetic noise radiating from such equi pment must 



ing th e set-up and wires (shielded twisted cable pairs) (jTardif et al. 



je suppressed by shield- 



2011 



Schoemaker et al. 



20081 ) . Moreover it is essential to have a rigid framewq r 



c. A mechanical resonance can occur 



in the cell/transducer system (at 70 Hz in 



Pengra et al. 



( 119991 )). and the noise ass ociated with 



mech anical vibration can be suppressed puting an additional mass to the frame ( jTardif et al. 



201ll). 



Once the oscillatory pressure is applied, the pressure must be measured. Most of the 
setups include piezoele ctric transduce r s to m easure the pressure difference over the capillary 
or the porous sample. 



Reppert et al 



response from 1 to 20 kHz. 



Tardif et al. 



(2001) proposed to use hydrophones that have a flat 



(1201 ll ) proposed to use dynamic transducers with a 



low- frequency limit 0.08 Hz and a maximum frequency of 170 kHz. 

The electrodes are u sually Ag/AgCl or platinium electrodes. The electrodes used by 



Schoemaker et al.l (120081 ) were sintered plates of Monel (composed of nickel and copper). 



The electrical signal must be measured using pre-amplifiers or a high-input impedance ac- 
quisition system. Since the impedance of the sample depends on the frequency, one must 
correct the measurement s from this varying- impedance to be able to have a correct stream- 



ing potential coefficient ( iReppert et al. 



200ll ) . Moreover the electrodes at top and bottom of 



the sample can behave as a c apacitor, requiring a correction using impedance measurements 



too (ISchoemaker et al. 



2008|). 



The sample is usually saturated and it is emphasized that the sample should be left until 
equilibrium with water. This equilibrium can be obtained by leaving the sample in contact 
with water for some time, and by flowing the water within the sample se yeral times by 



checking the pH and the water conductivity until an equilibrium is reached ( iGuichet et al. 



12 L. Jouniaux and C. Bordes 



20031 ). The procedure including water flow is better because the properties of the water can 
be measured. When the properties of the water are measured only before saturating the 
sample, the resulting water once in contact with the sample is not known. Usually the water 
is more conductive when in contact with the sample, and the pH can change. Recalling that 
the streaming potential is proportional to the zeta potential (which depends on pH) and 
inversely proportional to the water conductivity (eqJTJ), it is essential to know properly the 
pH and the water conductivity. 



4 MEASUREMENTS AND CALCULATIONS OF THE DYNAMIC 
ELECTROKINETIC COEFFICIENT 

The absolute magnitude of the st r eamin g potential coefficient normalized by the steady-state 



value was calculated by 



Packard! (119531 ) as: 



f(Xa) 



-iuit 



(20) 



which is equal to eq. [2j but expressed as a function of the parameter Y a = a^J^-, the 
transition frequency being obtained for Y a = 1 (Fig. [7]). The streaming potential coeffi- 
cient is constant up to the transition angular frequency, and then decreases with increasing 
fre quency. 



Sears fc Groves! (119781 ) measured the streaming potential coefficient on a capillary of 
radius 508 /im which was coated with clay-Adams Siliclad and then incubated with 1% 
bovine serum albumin, and filled with 0.02 M Tris-HCl at pH 7.32. They reported the 
streaming potential and the pressure difference as a function of frequency in the range — 20 
Hz. We calculated the resulting streaming potential coefficient (see Fig. [8]) which decreases 
from about 1.3x 10~ 7 to 4x 10~ 8 V/Pa. These authors computed the zeta potential and 
concluded that the zeta potential is independent of the frequency with an average value 
of 28.8 mV. Moreover they concluded that the zeta potential is also independent of the 
capillary radius and capillary length. 



The value of t 



Tardif et al. 
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le streaming potential coefficient on Ottawa sand measured at 5 Hz by 



(120111 ) was — 5.2x 10 7 V/Pa using a 0.001 mol/L NaCl solution to saturate the 



sample. Val ues between 1 and 2 x10 8 V/Pa were measured on samples saturated by 0.1 M/L 



NaCl brine (IPengra et al. 



19991 ) . A compilation of numerous streaming potential coefficients 



20101) 



measured on sands and sandstones at various salinities in DC domain ( lAUegre et al 
showed that C s o = — 1.2 x 10~ 8 aJ , where C s q is in V/Pa and o~f in S/m. A zeta potential of 
— 17mV can be inferred from these collected data, assuming the other parameters (see eq. d]) 
independent of water conductivity. These assumptions are not exact, but the value of zeta 
is needed for numerous modellings which usually assume the other parameters independent 
of the fluid conductivity. Therefore an average value of —17 mV for such modellings can be 
rat her exact, at least for medium with no clay nor calcite. 



Reppert et al.l (120011 ) calculated the real part and the imaginary part of the theoretical 
Packard's streaming potential coefficient (eq. [2]) for different capillary radii, (see Fig. [9]). It 
can be seen that the larger the radius is, the lower the transition frequency is, as shown above 
by the different theories. Recent developments by the group of Glover have been performed 
to build a new setup and to make further measurements on porous samples: two papers 



detail these studies in this special issue on Electrokinetics in Earth Sciences. 



5 CONCLUSION 

Since the theory of Pride in 1994, the dynamic behavior of the streaming potential is known 
for porous media. However few experimental results are avalaible, because of the difficulty to 
perform correct measurements at high frequency. Up to now, measurements of the frequency- 
dependence of the streaming potential have been performed up to 200 Hz on high-permeable 
samples. The main difficulty arises from electrical noise induced by mechanical vibration. 
Moreover it has been emphasized that the measurements must be corrected by impedance 
measurements as a function of frequency too because the impeda nce of the sample de 



pends on frequency. Further theoretical developments performed by 



Garambois fc Dietrich 



(120011 ) studied the low frequency assumption valid at frequencies lower than the transition 
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frequency. We show that this transition frequency, on a various collection of samples for 
which both formation factor and permeability are measured, is predicted to depend on the 
permeability as inversely proportional to the permeability. 
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Figure 1. The transition frequency f c = oj c /2tt (in Hz) predicted in the present study with ui c 
from eq. [12] with n = 10~ 3 Pa.s and Pf = W 3 kg/m 3 as a function of the permeability (in m 2 ). The 
transition frequency varies as logio(f c ) = —0.78logio(k) — 5.5. The parameters of the samples, F 
and fcg are measured from different authors on various samples cited in Tables 1, 2 and 3 




Figure 2. The sylphon bellows is driv en by a geophq ne-type push-pull driver to apply a sinusoidal 
motion to the sample, (modified from Packard! (|1953l )) 
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Figure 3. Experimental setup used by 



Thurston! (11952a! ) (modified from lThurstonl (|1952al )) 
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[figure 4. Experimental setup used by iGroves fc Sears! (|1975l ) (modified from 

ta3))- 



Groves &; Sears 
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Figure 5. Experimental setup used by 



osmosis measurements (modified from 



Pen gra et al 



1999) for streaming potential and electro- 



1999)) 
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Figure 6. Experimental setup used bv lReppert et al.l (|200ll ) (modified from 



Reppert et al 



(|200l|)). 
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Figure 7. The absolute magnitude of the normalized streaming potential coeffici ent calcu 



Packard! (j 19531 ) using eq.l20l where Y a = a y~^: equivalent to eq. [2] (modified from 



Packard 



atedjpy 
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Figure 8. The str eaming potential coefficient measured as a function of frequency by 
Sears Groves! (119781 ) on a capillary coated with clay, incubated with BSA in 0.02 M Tris-HCl. 
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. figu re 9. The real and imaginary part of the Packard's model (eqj2]) calculated bv lReppert et al. 
(|200ll ) for three capillary radii: 100/im(continuous line), 50;um(dashed line), 10/im(point line) (mod- 
ified from lReppert et alj ()200ll )). 



